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Abstract. We present the latest determination of the X-ray (2-10 keV) AGN luminosity 
function accounting for the selection effect due to X-ray absorption. The main results are: 
1) the inclusion of obscured AGN confirms the AGN differential luminosity evolution, but 
makes it less extreme than what is found selecting unobscured AGN in soft X-rays, and 
more similar to a pure luminosity evolution; 2) significant correlations are found between 
the fraction of obscured sources, the luminosity and the redshift, this fraction increasing 
toward both low AGN luminosities and high redshifts. We discuss our findings in a scenario 
for the formation and evolution of the structure in the Universe where the bulk of nuclear 
activity is produced at z~ 1 - 2. At the same redshifts also the star-formation rate reaches a 
maximum, and this age can therefore be regarded as the "golden age" for nuclear and galaxy 
activity. We discuss the current observational limits of this program and the improvements 
needed to obtain an unbiased census of the AGN and super-massive black hole (SMBH) 
population. 



1. Introduction 

Active Galactic Nuclei are not only witnesses 
of the phases of galaxy formation and/or as- 
sembly, but most likely among the leading ac- 
tors. Indeed, two seminal discoveries indicate 
tight links and feedbacks between SMBH, nu- 
clear activity and galaxy evolution. The first is 
the discovery of SMBH at the center of most 
nearby bulge dominated galaxies, and, in par- 
ticular, the steep and tight correlation between 
their masses and galaxy bulge properties (see 
e.g. Ferrarese & Merrit 2000, Gebhardt et al. 
2000, Marconi & Hunt 2003). The second dis- 
covery was originally due to the first deep X- 
ray surveys performed by ROS AT at the begin- 
ning of the 90' . They showed that the evolution 
of AGN is luminosity dependent. On average, 
the activity of Seyfert like objects rises up to za: 
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I and then decreases, while QSO activity rises 
smoothly up to z=2-3 or even further (Miyaji 
et al. 2000 and references therein). The former 
recalls the evolution of star-forming galaxies, 
while the latter recals the evolution of massive 
spheroids (Franceschini et al. 1999). However, 
soft X-ray surveys are biased against obscured 
sources, which, on the other hand, are very 
common in the local Universe. Indeed, the first 
imaging surveys above 2 keV obtained with 
ASCA and BeppoS AX found the first size-able 
samples of highly obscured AGN at z> 0.1 
and confirm, at least qualitatively, the predic- 
tions of standard AGN synthesis models for 
the Cosmic X-ray Background, CXB (Fiore 
et al. 1999, Akiyama et al. 2000, La Franca 
et al. 2002, Ueda et al. 2003). Chandra and 
XMM-Newton surveys confirm and expand 
this picture. On one side they resolved nearly 
100% of the CXB below 2 keV (Giacconi et 
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al. 2002, Brandt et al. 2001) and confirm the 
strong luminosity dependent density evolution 
of soft X-ray sources (Hasinger 2003, 2005). 
On the other side deep and large area surveys 
up to 10 keV clearly showed that AGN activ- 
ity spans a range of optical to near-infrared 
properties much greater than it was thought 
based on optically and soft X-ray selected 
AGNs (Fiore et al. 2003, Koekemoer et al. 
2004, Barger et al. 2005). The main reason 
is that hard X-ray selection provides a more 
complete and direct view of AGN activity, be- 
ing less biased than optical or soft X-ray se- 
lection against obscured sources. For exam- 
ple, 2-10 keV surveys pick up AGN relatively 
bright in X-rays but with extremely faint opti- 
cal counterparts. The majority of these sources 
have been identified as highly obscured, high 
luminosity AGN at z; 1, the so called type 2 
QSOs (~ 20 in the HELLAS2XMM survey, 
3/4 confirmed through optical spectroscopy, 
Fiore et al. 2003, MignoU et al. 2004, Maiolino 
et al. 2006, Cocchia et al. 2006, about 30 in 
the CDFSh-CDFN, half of which confirmed 
through optical spectroscopy; a dozen in the 
CLASXS survey, Barger et al. 2005). At the 
opposite of the X-ray to optical flux ratio dis- 
tribution, Chandra and XMM-Newton hard X- 
ray surveys discovered moderately obscured 
sources with AGN luminosity, in otherwise 
inactive, optically bright, early type galaxies 
(named X-ray bright, optically normal galax- 
ies, XBONGs, Fiore et al. 2000, Comastri et 
al. 2002). 

In this paper we discuss how the inclusion 
of obscured AGN affects the determination of 
the AGN luminosity function, and discuss our 
findings in the framework of semi-analytical 
models for the formation and evolution of the 
structure in the Universe. 

2. The evolution of hard X-ray 
selected sources 

Figure [l] gives an overview of the flux limits 
and surveyed areas of major AGN surveys car- 
ried out over the last years in the 2-10 keV 
band. In this paper we use the source samples 
given in Table 1 , which include the deepest sur- 
veys performed with Chandra as well as larger 
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Fig. 1. Solid angles and flux limits of AGN surveys 
carried out in the 2-10 keV band. Triangles repre- 
sent serendipitous surveys constructed from a col- 
lection of pointed observations (Chandra open sym- 
bols; XMM-Newton filled symbols; BeppoSAX and 
ASCA surveys are also reported). The asterisks rep- 
resent surveys covering contiguous areas. 



area surveys performed with XMM-Newton 
and ASCA, plus the so-called Piccinotti sam- 
ple of local AGN. 

We estimated the 2-10 keV luminos- 
ity function by fitting the expected number 
of AGN in bins of luminosity, redshift and 
rest frame absorbing column density Nh (La 
Franca et al. 2005). This allows us to take into 
account observational selection effects. In par- 
ticular, we correct for the selection effect due to 
X-ray absorption, and for the incompleteness 
of the optical spectroscopy identification (see 
La Franca et al. 2005 for details). Our results 
extend those of Fiore et al. (2003), Cowie et al. 
2003 and Barger et al. (2005). In all these three 
papers no correction for the X-ray absorption 
is adopted. Fiore et al. (2003) assign a redshift 
to the sources without spectroscopic identifi- 
cation using statistical arguments; Cowie et al. 
(2003) and Barger et al.(2005) estimate upper 
limits to the AGN density by assigning to the 
unidentified sources the redshifts correspond- 
ing to the centers of each Lx-z bin. 

Figure |2 and |3 show the best fit Nh 
distribution and best fit luminosity functions 
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Table 1. 2-10 keV surveys 



Sample 


Tot. Area. 


Flux limit 

10"'^ cgs 


# sour. 


% z-spec 


HELLAS2XMM 


1.4 


6.0 


232 


70% 


rnFN faint" 


0.0369 


1.0 


95 


59% 


CDFN bright* 


0.0504 


3.0 


51 


65% 


CDFS faint" 


0.0369 


1.0 


75 


62% 


CDFS bright* 


0.0504 


3.0 


52 


60% 


Lockman Hole^ 


0.126 


2.6 


55 


75% 


HBS28 


9.8 


22 


28 


100% 


AMSSn 


69 


30 


74 


100% 


Total 






662 


75% 



" Inner 6.5 arcmin radius; * outer 6.5-10 arcmin annulus; inner 12 arcmin radius; inner 4.5 arcmin radius. 
See La Franca et al. 2005 for details. 



for a luminosity-dependent density evolution 
(LDDE) model (La Franca et al. 2005, see also 
Miyaji et al. 2000 for a similar parameteriza- 
tion). The parameters of the evolving luminos- 
ity function and of the dependencies of the Nh 
distribution by L(2-10 keV) and z have been 
fitted simultaneously. The intrinsic Nh distri- 
bution (dotted lines in figure |2j is flat above 
10^' cm"^, while the fraction of objects with 
Nh < 10^' cm"^ is one of the model param- 
eters. The dashed lines show the predictions 
when the selection effects due to X-ray absorp- 
tion, which pushes sources below the X-ray de- 
tection limit, and to the incompleteness of the 
spectroscopic identification are taken into ac- 
count. As expected, most AGN with Nh > lO^'' 
cm"^ are lost in even the deepest Chandra and 
XMM 2-10 keV surveys. In addition, note that 
about 2/3 of the AGN with lO^"* <Nh < lO^^ 
cm"^ are also lost at z< 1. Our best fit lumi- 
nosity function recovers all these highly X-ray 
obscured AGN, as well as optically obscured 
AGN, whose optical counterpart is too faint to 
allow a spectroscopic identification. 

The best fit 2-10 keV luminosity function 
is compared to the best fit 0.5-2 keV total lu- 
minosity function of Miyaji et al. (2000) and 
to flie best fit 0.5-2 keV type 1 AGN luminos- 
ity function of Hasinger et al. (2005) in figure 
|3] These functions falls shorter than the 2-10 
keV luminosity function by a factor 3-10 at 
luminosities logL^ = 42.5 - 43 in the three 
lowest redshift bins. (At z> 1.5 the compari- 
son is less informative because the present 2- 
10 keV data do not go deep enough to provide 
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Fig. 2. N/y distributions in four luminosity and red- 
shift bins. The dotted lines are the assumed N/y dis- 
tributions, the dashed lines are the expectations tak- 
ing into account the selection effects and the contin- 
uous lines are the observed distributions. 



samples of low luminosity AGNs large enough 
to constrain adequately their space density.) 
Conversely, a better agreement between the 
0.5-2 keV and the 2-10 keV luminosity func- 
tions is found at the highest luminosities sam- 
pled, with the exception of the z=l-1.5 bin. 
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Fig. 3. The 2-10 keV luminosity function in four redshift bin, compared to the 0.5-2 keV luminosity func- 
tion of Miyaji et al (2000, left panel, dotted lines), with the 0.5-2 keV type 1 AGN luminosity function of 
Hasinger et al. (2005, left panel, solid lines) and with the optical (B band) luminosity function of Croom 
et al. (2004, right panel). Two optical luminosity function are plotted in each quadrant, corresponding to 
redshift bridging the range used for the 2-10 keV luminosity function. The conversion factor to pass from 
the 0.5-2 keV to the 2-10 keV band has been calculated assuming a power law spectrum with a = 0.8 while 
that to pass from the B band to the 2-10 keV band has been computed following Marconi et al. 2004. 



Hasinger et al. (2005) exclude from their anal- 
ysis X-ray (and/or) optically obscured AGN, 
which in any case are not an important popula- 
tion in 0.5-2 keV surveys. Indeed the Hasinger 
et al. luminosity function is nearly identical to 
the Miyaji et al. (2000) one at least at z< 1.5. 
(At z> 1.5 the Miyaji et al. 2000 best fit over- 
predicts the number of low luminosity AGN, 
probably because to the larger uncertainties at 
high redshift and low luminosity due to the 
much shallower data used in comparison with 
Hasinger et al. 2005). Obscured AGNs are re- 
covered in the 2-10 keV luminosity function. 
The result is that the evolution of this lumi- 
nosity function deviates less from a pure lumi- 
nosity evolution than the soft X-ray luminosity 
functions. At least part of the strong luminos- 
ity dependent density evolution claimed based 
on the 0.5-2 keV data is therefore due to the 
exclusion of obscured AGN. 

Figure |3l compares the 2-10 keV and op- 
tical luminosity functions. Luminosity depen- 
dent conversion factors to pass from the B band 
to the 2-10 keV band have been computed 
following Marconi et al. (2004). It should be 



borne in mind that the uncertainty on the con- 
version factor may be as large as a factor of 
two. Note that the low luminosity end of the 
optical luminosity function is always 1, 1.5 
dex higher than that of the X-ray luminosity 
function. This is due to the difliculty in se- 
lecting low luminosity AGNs against their host 
galaxy in the optical band. The optical lumi- 
nosity function is similar or slightly lower than 
the 0.5-2 keV luminosity function at z> 0.5, 
and therefore the same comments given above 
for the soft X-ray luminosity function apply 
also in this case: e.g. a large fraction of ob- 
scured accretion may be missed by both optical 
and soft X-ray selection. 

Figure |3 shows the fraction of X-ray ob- 
scured AGN (Nh > 10^^ cm"^) as a function 
of the 2-10 keV luminosity and the redshift. 
The long dashed lines are the best fit intrin- 
sic distributions while the short dashed lines 
are the expectations taking into account all se- 
lection eff'ects described above. Both the ob- 
served and best fit fraction of obscured AGN 
at z<l decrease strongly with the AGN lumi- 
nosity, a behavior already noticed in the lit- 
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erature since the first Einstein systematic ob- 
servations of QSO (Lawrence & Elvis 1982), 
and confirmed quantitatively by Ueda et al. 
(2003). Note that this trend becomes more ev- 
ident putting together deep surveys, well sam- 
pling low luminosity AGN, and large area sur- 
veys, which can provide large samples of lumi- 
nous AGN. Note also that selection effects af- 
fects in a similar way low and high luminosity 
AGN: the ratio of the intrinsic and predicted 
fractions of obscured AGN is nearly constant 
with the luminosity. This is not the case when 
considering the dependence of the fraction of 
obscured AGN with the redshift. The reason 
is that obscured AGN are much more likely 
to be missed in X-ray surveys at low redshift 
than at high redshift, because the photoelectric 
cut-off is quickly redshifted toward low X-ray 
energies. Indeed, comparing the best fit intrin- 
sic distribution to the distribution expected af- 
ter the selection effects, we note that while we 
loose just about 15% of the obscured AGN at 
T-2 slightly less than half are lost at z< 0.2. By 
correcting for this selection effect we find that 
the intrinsic fraction of obscured AGN in the 
luminosity range lO'*-' - 10"^^ still increases by 
a factor 1 .8 from z=0 to z=2. This new result is 
due a better coverage of the luminosity -redshift 
diagram in comparison with previous work. 

Figure |5] shows the luminosity-redshift 
plane for the surveys in Table 1. The arrows 
indicate the directions in which the fraction of 
obscured AGN increases, according to figure 
m Excluding the CDFS and CDFN data would 
strongly limit the number of AGN with logL(2- 
10)< 44 (those more likely to be obscured, ac- 
cording to the left panel of figure 0}, making 
difficult to asses any trend with the redshift. 
The arrows indicate also the direction in which 
the flux decreases. Indeed a strong correlation 
of the number of obscured AGN and the flux 
has been reported in the past (Piconcelli et al. 
2003, Ueda et al. 2003, Perola et al. 2004) and 
it is expected by AGN synthesis models of the 
CXB (Comastri et al. 2001). Figure |5] shows 
this correlation for the source sample in Table 
1, along with the best fit LDDE model of La 
Franca et al. (2005). In conclusion, the corre- 
lations found between the fraction of obscured 
AGN and luminosity and redshift of figure |4] 




42 43 44 45 46 0.5 1 1.5 2 2.5 3 3.5 4 
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Fig. 4. Observed fraction of X-ray obscured (N/^ > 
10^- cmr^) AGN as a function of L(2-10keV) and z. 
The long dashed lines are the best fit intrinsic distri- 
butions. The short dashed lines are the expectations 
taking into account all selection effects. 



confirm and extend previous determinations, 
based on smaller and shallower surveys. 

3. Discussion 

Our determination of the 2-10 keV AGN lu- 
minosity function, accounting for selection ef- 
fects due to nuclear obscuration by gas and 
dust, confirms the AGN differential luminosity 
evolution, but makes it less extreme than what 
is found selecting unobscured AGN only (see 
e.g. Hasinger et al. 2005). This is important 
for both models that make use of the AGN lu- 
minosity function (to reproduce the X-ray and 
IR Cosmic backgrounds for example), and for 
models which try to explain the AGN lumi- 
nosity function, as the semi-analytic, hierarchi- 
cal clustering model proposed by Menci et al. 
(2004). Figure |6l compares the AGN number 
density as a function of z with the predictions 
of the Menci model. The result is qualitatively 
similar to that reported by Menci et al. (2004). 
The trend of lower luminosity AGN peak- 
ing and increasingly lower redshift is found 
in both data and model, which however pre- 
dicts a number of low-to-intermediate lumi- 
nosity (Seyfert like) AGN at z= 1.5-2.5 higher 
than what is observed by a factor of a few. This 
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Fig. 5. Left panel: the luminosity-redshift plane for the surveys in Table 1. Crosses indicate AGN with 
Nh > 10"^ cm"^. The arrows indicate the directions in which the fraction of obscured AGN increases, 
according to figure |4] Right panel: observed fraction of X-ray obscured (N/y > 10^" cm"-) AGN as a 
function of F(2-10keV) 



disagreement can be due to at least two broad 
reasons: either La Franca et al. (2005) under- 
estimate the number of highly obscured AGN 
missed at z= 1.5-2. 5 in Chandra and XMM- 
Newton surveys, or there are problems with the 
model. 

About the first possibility, it must be noted 
that most obscured AGN selected below 10 
keV have column densities in the range Nh ~ 
]^q22-23 ^jjjj Qjjjy ^ handful of the faintest 

sources which may be Compton thick (Nh > 
10^4 cm-^, just ~ 4% in the CDFS, Tozzi 
et al. 2006). So we still may be viewing 
just the tip of the iceberg of highly obscured 
sources. Our estimates of the number of ob- 
scured AGN missed in today X-ray surveys 
are based on large extrapolations from what 
we know about the fraction of obscured AGN 
in the local Universe. Compton thick objects 
may well be more common at high redshift, 
as suggested, for example, by Fabian (1999) 
and Silk & Rees (1998). An alternative ap- 
proach to find Compton thick AGN at z>l 
is to select sources with AGN luminosities 
in the mid-infrared and faint near-infrared 
and optical emission (Martinez-Sansigre et al. 
2005). These authors estimate that probably 
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Fig. 6. The evolution of the number density of 
AGNs selected in the 2-10 keV band in three bins of 
luminosity (43 < logLx < 44 44 < logLx < 44.5 
44.5 < logLx compared to the prediction of the 
model of Menci et al. (2004) 



more than half of high luminosity QSOs are 
highly obscured, although with quite large un- 
certainties. Unfortunately, the X-ray proper- 
ties of these infrared selected sources are not 
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Fig. 7. A working scenario for Compton thin AGN 



known, and therefore it is difficult to under- 
stand how the mid-infrared selection compares 
with the X-ray one. In particular, it is not clear 
which is the fraction of the mid-infrared se- 
lected type 2 AGN which would have been 
selected by X-ray surveys. Answers to these 
questions will soon come from the study of 
fields with both X-ray and mid-infrared cov- 
erage (e.g. the ELAIS-Sl field, Puccetti et al. 
2006, Feruglio et al. 2006 in preparation, and 
the COSMOS fields), and from deep X-ray 
follow-up observations of the mid-infrared se- 
lected sources in the Spitzer First Look Survey. 

To avoid any possible selection effect, for 
an unbiased census of the AGN population 
making the bulk of the CXB and an unbi- 
ased measure of the AGN luminosity func- 
tion at z=l-2, sensitive observations extend- 
ing at the peak of the CXB are clearly needed. 
More specifically to resolve ~ 50% of the 
CXB in the 20^0 keV band we need to go 



down to fluxes of 10"'^ erg cm"^ s"' in this 
band. This can be achieved only by imaging 
X-ray telescopes, (see e.g. Fiore et al. 2004, 
Ferrando et al. 2005). Key issues are: a) high 
collecting area; b) sharp PSF (15 arcsec or less 
Half Energy Diameter); c) low detector inter- 
nal background. 

About models, as in all complex processes, 
there are several areas which may be critical. 
For example, the prescriptions adopted by the 
Menci model to switch an AGN on and to com- 
pute its feedback on the host galaxy may be 
too simple, or, more in general, the descrip- 
tions of the mechanisms regulating the amount 
of cool gas in low-mass host galaxies and the 
physical mechanism at work at small accre- 
tion rates may be inadequate, as well as the 
statistics of DM condensations. Further con- 
straints to this model, which may shed light on 
these issues, come from the observed correla- 
tions between the fraction of obscured AGN 
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with luminosity and redshift. These correla- 
tions are at odds with popular AGN unified 
Schemes (see e.g. Lamastra et al. 2006), and 
may suggest that low luminosity, Seyfert like 
AGN and powerful QSOs are intrinsically dif- 
ferent populations, with different obscuration 
properties, caused by different formation histo- 
ries, a bimodal behaviors reminding that of the 
color distribution of galaxies (see e.g. Menci 
et al. 2005). The Seyfer-like object popula- 
tion could be due to nuclear activation during 
loose galaxy encounters (fly-by, Cavaliere & 
Vittorini 2000) at sub Eddington levels, and the 
second population could be due to nuclear ac- 
tivation during major mergers, in the process 
of galaxy assembly. Seyfert like AGN could 
be mostly associated to galaxies with merging 
histories characterized by small mass progen- 
itors while QSOs may be associated to large 
mass progenitors, as sketched in the cartoon 
of figure In the first case nuclear accretion 
and star-formation could be self-regulated by 
feedbacks, which can therefore be effective in 
leaving available cold gas that can both cause 
an obscuration of the nucleus and be accreted 
during subsequent galaxy encounters. In these 
galaxies gas and dust lanes can efficiently ob- 
scure the nucleus along many lines of sight (a 
scenario similar to that outlined by Matt 2000). 
In the second case case feedback could be less 
effective in reheating/expelling the cold gas, 
most of which is rapidly converted in stars at 
high z. The obscuration properties of the two 
populations could be different in terms of gas 
geometry, covering factor, density, ionization 
state, metallicity, dust content and composi- 
tion. A quantitative comparison between the 
prediction of the Menci model about the frac- 
tion of obscured AGNs with the observation is 
in progress. This comparison can help in both 
understanding which is the leading physical 
mechanism responsible for the activation and 
obscuration of AGNs of different luminosities, 
and in understanding the role of relative feed- 
backs between nuclear activity, star formation 
and galaxy evolution, as a function of the host 
mass and luminosity. 
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